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The results on directed and elliptic flow for Pb+Pb at the full energy of the SPS (158 GeV/A) and 
from the first year of Au + Au at RHIC (^/s NN = 130 GeV) are reviewed. The different experiments 
agree well and a consistent picture has emerged indicating early time thermalization at RHIC. 



o 
o 

(N 



> 

o 
o 



o 

■ft 

=5 



X 



Elliptic flow from non-central collisions has become a 
prime hadronic signature for pressure at early times in 
relativistic heavy ion collisions jl], |J. Parton Cascade 
Model calculations |!| show that the elliptic flow builds 
up in the first few fm/c and then remains constant. In 
this model the amount of elliptic flow produced is pro- 
portional to the assumed parton-parton scattering cross 
section. Thus it is thought that rescattering converts 
the initial space anisotropy of the overlap region of non- 
central collision to the momentum anisotropy of elliptic 
flow. As the initial lens-shaped overlap region expands it 
becomes more spherical, quenching the driving force that 
produces the elliptic flow. Thus elliptic flow is sensitive 
to the number of interactions and is considered to be a 
measure of the degree of thermalization at early time. 



I. INTRODUCTION 

Directed and elliptic flow are defined by the v\ and 
V2 coefficients in the Fourier expansion of the azimuthal 
distribution of particles relative to the reaction plane: 

1 + 2vi cos((f> - * r ) + 2v 2 cos(2(0 - * r )), (1) 

where <fi is the azimuthal angle of a particle and \& r is the 
azimuthal angle of the reaction plane. The coefficients 
are usually evaluated differentially as a function of y, pt, 
and centrality. 

Three methods of analysis will be discussed. In the 
conventional method 0] of correlating particles with an 
event plane the flow coefficients are given by: 



produces the squares of the coefficients, so that one has 
to take the square root of the correlation effect: 



,,obs 



(cos(n((j>i - ^„))), 



(2) 



where ^ n is the observed event plane of order n. Since 
the observed event plane is not the true reaction plane, 
the observed coefficients have to be corrected by divid- 
ing by the resolution of the event plane, which is esti- 
mated by measuring the correlation of the event planes 
of subevents. 

The flow coefficients can also be obtained by the pair- 
wise correlation || of all the particles without refering 
to an event plane. This two-particle correlation method 



(3) 



Recently, a multiparticle correlation method has been 
proposed H using cumulants. For four particles, one 
calculates the four-particle correlation minus twice the 
square of the two-particle correlations: 



+ (cos(n(0i - <j> 3 ))) (cos(n(0 2 - M)) 
+ (cos(n(0i - (j) 4 ))) (cos(n(0 2 - <fe)))- 



(4) 



Here one must take the fourth root of the result, and 
naturally the statistical errors are larger than for the 
two-particle analysis. However, this method has the ad- 
vantage of eliminating two-particle non-flow effects which 
are not correlated with the reaction plane, such as HBT 
and resonance decay. 



II. SPS 
A. NA49 
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FIG. 1: NA49 directed and elliptic flow for pions from a min- 
imum bias trigger as a function of rapidity. The open points 
are reflected. The data are still preliminary. 
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From NA49, initial flow results from part of the data 
set were published 0, |) a few years ago. The full set of 
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data has now been analyzed Q. The Vi and v-i values 
for pions as a function of rapidity are shown in Fig. [l]. 
The values have been integrated over p t by taking mean 
values weighted with the cross section. The v\ values 
have been corrected for conservation of momentum flio| . 
Before this correction, v\ crossed mid-rapidity at +0.6%, 
while after the correction, with no adjustable parameters, 
it can be seen that the curve now crosses mid-rapidity at 
zero. Other non-flow effects [[ll], [l2| have been verified 
to be small using the following tests. Short-range cor- 
relations, like HBT, were determined to be not large by 
using striped sub-events where only physically separated 
particles were correlated. Coulomb effects were deter- 
mined to be not large by looking at tt + and 7r~ particles 
separately. Neutral resonances, like p°, were shown to 
be not important by correlating tt~ with an event plane 
determined only by other tt~ particles. The pion flow 
values integrated over both p t and y are show in Fig. ^ 
as a function of centrality bin. 



£ -5 
^-10 
-15 
-20 
-25 
-30 



i 


— 














• • • • • 






r pions 












i 


— 








>- 



2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 7 

rapidity 

FIG. 3: NA49 (circles) and WA98 (squares) pion data are 
compared. Plotted is the directed flow versus rapidity. The 
vertical lines are at mid-rapidity and beam rapidity. The data 
are preliminary. 
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FIG. 2: NA49 directed and elliptic flow for pions as a function 
of centrality bin. The centrality increases to the left. The data 
are still preliminary. 
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FIG. 4: NA49 (circles) and NA45 (squares) charged particle 
data are compared. Plotted is the elliptic flow at mid-rapidity 
versus the centrality. The centrality increases to the left be- 
cause the abscissa is the cumulative most central percentage 
of all events that are in the bin. The data are preliminary. 



WA98 has measured |l3| the directed flow of pions near 
target rapidity in the Plastic Ball. These data are plotted 
in Fig. U reflected into the forward hemisphere, together 
with the NA49 data, ft can be seen that at projectile ra- 
pidities the directed flow is huge. WA98 also has shown 
that for protons v\ is just as large, but in the positive di- 
rection. This indicates that looking at projectile rapidity 
fragments might be a good way to determine the reaction 
plane. However, even though v\ is large, the resolution 
of the event plane also depends on the square root of the 
multiplicity, which in this rapidity region is low. 



Using silicon drift chambers NA45 has measured |T^] 
the flow coefficients for charged particles. In Fig. |J their 
elliptic flow values at mid-rapidity are plotted versus cen- 
trality, together with the same quantity for pions from 
NA49. The agreement is good. NA45 also has a RICH 
detector which identifies high p t pions. When correlated 
with an event plane determined from their TPC, the pt 
dependence falls nicely on a line extrapolated from the 
lower p t results of NA49. However, when they correlate 
pairs of high p t pions from the RICH using Eq. ||, their 
results are somewhat higher p5[|. 
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D. NA50 



B. PHENIX 



NA50 has analyzed the neutral transverse energy from 
their calorimeter in the rapidity region from 1.1 to 2.3. 
NA50 has shown elliptic flow as a function of central- 
ity [jl6|, but it is hard to compare their results with 
the other experiments because of the different method 
of analysis. 



III. RHIC 
A. STAR 

STAR results for elliptic flow for charged particles JTtJ 
and identified particles have been published. The 
flow is larger than at the SPS and comes close to the 
hydrodynamic predictions ]l9| for the mid-central and 
central data. The identified particle flow shows a varia- 
tion with mass as predicted by hydrodynamics. However, 
the high p t behavior [|o) exhibits a saturation effect well 
below where hydro would saturate. It has been suggested 
that at high p t elliptic flow could be due to parton energy 
loss (jet quenching) pj] . 

The doubly integrated elliptic flow values as a function 
of centrality are shown in Fig. ||, using the conventional 
event plane method (Eq. ||). The asymmetric error bars 
indicate the systematic uncertainty due to possible non- 
flow effects. A four-particle correlation analysis |2^] using 
Eq. ^ is also show in the figure. Indeed, when two-particle 
non-flow effects are eliminated, the points do come down 
to about the limits of the systematic errors. It should be 
noted that part of the decrease could be due to fluctua- 



The PHENIX results (24) for charged particle elliptic 
flow as a function p t for a minimum bias trigger is shown 
together with the STAR results in Fig. |. PHENIX uses 
the pair-wise method of Eq. ^[ The agreement is ex- 
cellent. However, at the moment, the doubly-integrated 
elliptic flow values as a function of centrality are about 
50% higher for PHENIX than for STAR. Only about half 
of this discrepancy can be accounted for by the higher p t 
threshold of PHENIX. 




tion effects 23 which also would go in the same direction. 



However, the agreement with hydro for mid-central and 
central collisions is not affected. 



FIG. 6: STAR (circles) and PHENIX (squares) charged par- 
ticle data for a minimum bias trigger are compared. Plot- 
ted is the elliptic flow versus the transverse momentum. The 
PHENIX data are preliminary. 
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STAR Au + Au at \| s NN = 130 GeV: 
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FIG. 5: STAR elliptic flow for charged particles versus cen- 
trality. The results using the conventional event plane method 
(circles) and the four-particle cumulant method (squares) are 
shown. The later are preliminary. 



The PHOBOS elliptic flow values [£5| as a function 
of centrality are shown in Fig. [?], together with the 
STAR values. The agreement is excellent. Fig. || shows 
PHOBOS elliptic flow for a minimum bias trigger mea- 
sured over a wide range of pseudorapidity. Where there 
are STAR measurements, in a narrow range near mid- 
rapidity, the agreement is good with the flat STAR dis- 
tribution. A 3D hydro calculation |^(| gets rather flat 
distributions over a wide range of pseudorapidity for all 
reasonable initial conditions. Thus the fast fall off away 
from midrapidity in Fig. || is not understood at present 
in a model with thermalization. 



IV. CONCLUSIONS 

Good agreement between the various experiments has 
been achieved. There might be some indication that the 
pair-wise method for high p t particles indicates non-flow 
effects. The four-particle method at RHIC indicates that 
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FIG. 7: STAR (circles) and PHOBOS (squares) charged par- 
ticle data are compared. Plotted is the elliptic flow versus 
centrality. The PHOBOS data are preliminary and the error 
bars are only statistical. 



there are non-flow effects for peripheral collisions. The el- 
liptic flow increases from the SPS to RHIC, and at RHIC 
is close to the hydro limit. This is taken to be a signature 
for early time thermalization at RHIC. 



0.08 
0.07 
0.06 
0.05 
0.04 - 
0.03 — 
0.02 — 
0.01 



-6 



4+ 



-2 



♦++■ 



2 4 6 
Pseudorapidity 



FIG. 8: PHOBOS charged particle elliptic flow for a min- 
imum bias trigger plotted versus pseudorapidity. The data 
are preliminary and the error bars are only statistical. 
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